Stability and structure of atomic chains on Si(lll) 
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We study the stability and structure of self-assembled atomic chains on Si(lll) induced by mono- 
valent, divalent and trivalent adsorbates, using first-principles total-energy calculations and scanning 
tunneling microscopy. We find that only structures containing exclusively silicon honeycomb or sili- 
con Seiwatz chains are thermodynamically stable, while mixed configurations, with both honeycomb 
and Seiwatz chains, may be kinetically stable. The stability and structure of these atomic chains 
can be understood using a surprisingly simple electron-counting rule. 



I. INTRODUCTION 

Understanding and controlling the structure of sur- 
faces on the atomic level is of tremendous technologi- 
cal importance. This is especially true for the growth 
of semiconductor nanostructures, where the competition 
between thermodynamics and kinetics can play a deci- 
sive role. The reconstruction of semiconductor surfaces 
is driven by the elimination of dangling bonds and the 
minimization of surface stress, with a striking diversity 
of outcomes. Despite this variety, even very elaborate 
architectures are generally comprised of a small num- 
ber of elementary structural building blocks. Dimcrs 
and adatoms on Si(100) and Si(lll), respectively, are 
the best known strategies for reducing the number of 
dangling bonds. Tetramers and pentamers, encountered 
on Si(114)i, Si(113)£, Si(110)M, and Si(331)£ constitute 
more complex units. For adsorbate-induced surface re- 
constructions of the Si(lll) surface, honeycomb^^ and 
Seiwatz chains 9 have recently emerged as universal build- 
ing blocks^ (see Fig. [T|). These form the basis of a large 
class of atomic chain reconstructions which have been 
the focus of intense research because of their fascinat- 
ing quasi one-dimensional electronic propertie a 11 ' 12 ! 13 ' 14 . 
The fact that only silicon atoms participate in the for- 
mation of the honeycomb and Seiwatz chains means that 
a number of different adsorbates can induce a chain re- 
construction, simply by donating the correct number of 
electrons to the substrate. 

We have recently revealed remarkable systematics in 
this class of adsorbate-induced reconstructions, relat- 
ing the valence state of the adsorbate to the allowed 
coverages and periodicities of the resulting adsorbate 
chains—. All experimentally observed phases satisfy a 
simple electron-counting rule: the adsorbates must pro- 
vide either one electron (to stabilize a (3x1) honeycomb- 
chain unit), or two electrons (to stabilize a (2x 1) Seiwatz- 
chain unit)»i£ For monovalent adsorbates the only exper- 
imentally observed configuration is obtained with an ad- 
sorbate coverage of 1/3 MLj^ This corresponds to one 
adsorbate donating one electron per (3x1) honeycomb 
chain unit as shown in Fig. [TJa) , in agreement with the 



FIG. 1: (Color online) Top (upper) and side (bottom) view 
of the two prototypical silicon chains separated by channels 
that accommodate the adsorbate atoms (filled blue circles), 
(a) Honeycomb chain with (3x1) unit cell, (b) Seiwatz chain 
with (2x1) unit cell, (c) Mixed chain structure with (5x1) 
unit cell. The (lxl) unit cell of the (lll)-oriented substrate 
is indicated in grey. The crystallographic directions of the 
substrate are also indicated. 



electron-counting rule. For divalent adsorbates, only 1/6 
ML (half the monovalent adsorbate coverage) is required 
to stabilize honeycomb chains, since a divalent adsorbate 
donates two electrons. The pure Seiwatz chain structure 
shown in Fig. |Tfb) is stable for 1/2 ML of divalent ad- 
sorbates, resulting from two donated electrons per (2x1) 
Seiwatz chain unit. A mixed chain phase, alternating 
between honeycomb and Seiwatz chains as shown in Fig. 
[He), observed for divalent adsorbates, at intermediate 
coverage of 3/10 ML, also satisfies the electron-counting 
rule. Here the divalent adsorbates supply one electron 
to the honeycomb chain unit and two electrons to the 
Seiwatz chain unit. Similarly, trivalent adsorbates at a 
coverage of 2/10 ML donate the same number of elec- 
trons, allowing the same mixed phase to be stabilized. 

In this report we subject these observed systematics, 
and the electron-counting rule deduced from them, to 
more detailed theoretical scrutiny. Specifically, we ex- 
amine theoretically the thermodynamic and kinetic sta- 
bility of several chain reconstructions of Si(lll) using 
first-principles total-energy methods. We compare the 
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resulting pictures that emerge for three prototypical ad- 
sorbates: monovalent (Na), divalent (Ca), and trivalent 
(Gd). For each adsorbate we compare a number of re- 
constructions with different adsorbate coverage, includ- 
ing the bare Si(lll)-(7x7) reconstruction and, for Ca 
and Gd adsorbates, the silicide phase experimentally ob- 
served at higher coverage. For each candidate reconstruc- 
tion we determine the surface energy as a function of 
adsorbate chemical potential, and thereby determine the 
energy ordering of our candidate reconstructions at any 
thermodynamically allowed value of chemical potential. 
We find that the thermodynamically stable chain recon- 
structions are formed exclusively from either honeycomb 
or Seiwatz chains, with mixed phases slightly higher in 
energy. We also argue that for Gd adsorbates the exper- 
imentally observed mixed phase, which combines honey- 
comb and Seiwatz chains, while not thermodynamically 
stable, is kinctically stable. 



II. METHODS 

We used first-principles total-energy calculations to de- 
termine equilibrium geometries and relative surface ener- 
gies. The calculations were performed in a slab geometry 
with up to six layers of Si plus the reconstructed surface 
layer. All atomic positions were relaxed except for the 
bottom layer, which was passivated. Total energies and 
forces were calculated within the generalized-gradient ap- 
proximation to density- functional theory (DFT), using 
projector-augmented- wave (PAW) potentials 17 ' 18 . We 
checked that the slab thickness, plane-wave cutoff, and 
sampling of the surface Brillouin zone were each suffi- 
cient to converge the relative surface energies to within 
1 meV/A 2 . 

For calculations with Gd adsorbates, the seven 4/ elec- 
trons were treated explicitly as valence states. The possi- 
bility of magnetic order among Gd atoms within a single 
fully occupied channel was investigated in one case, with 
the result that ferromagnetic ordering was slightly pre- 
ferred, by 0.1 eV, to antiferromagnetic ordering. Based 
on this finding we assumed ferromagnetic order for all Gd 
phases. We also found that putting the 4/ electrons in 
the core led to only insignificant changes to the calculated 
absolute surface energies. This establishes that magnetic 
order among the Gd atoms plays no substantive role in 
the stability of different surface phases. 

Growth and STM experiments were carried out in a 
ultra-high vacuum chamber with a residual gas pressure 
of 3 x 10~ n mbar equipped with an Omicron LT-STM. 
Boron-doped Si(lll) with a resistivity of 5 fi-cm was 
heated by passing a direct current. Gd was evaporated 
from a water-cooled e-beam evaporator. 

We consider first the stability of atomic chains in- 
duced by monovalent adsorbates such as the alkali met- 
als Li, Na, K, Rb, and Cs. These are known to induce a 
chain reconstruction on Si(lll) exhibiting simple (3x1) 
periodicity. The widely accepted structural model, the 



so-called honeycomb chain-channel (HCC) model-i^, is 
shown in Fig. [lja). The HCC model consists of Si honey- 
comb chains aligned along the [110] direction, separated 
by empty channels. The adsorbate atoms occupy sites 
within these channels. 

The presence of an adsorbate, such as Na, in the HCC 
model does not allow direct comparison between its sur- 
face energy and the surface energy of clean reconstructed 
Si(lll)-(7x7). The proper way to compare energies of 
structures differing in stoichiometry is via the chemical 
potentials /isi and /iNa of the constituents^, which are 
the energy per atom available in the reservoirs with which 
the surface is assumed to be in equilibrium. The surface 
energy (per unit area) is then 

7 = Esaxf/A = (E tot - nsi/^Si - n Na ^ Na )/A, (1) 

where Etot is the total energy of a double-sided slab 
whose unit cell, with total area A, contains nsi Si atoms 
and riNa Na atoms. Since the surface is in equilibrium 
with the bulk Si substrate, nsi is the energy per atom in 
bulk Si. The adsorbate chemical potential, /iNa, however, 
corresponds to a real physical variable that can be exter- 
nally tuned by, for example, varying the partial pressure 
of Na. Intuitively, increasing the Na partial pressure will 
increase the stability of structures with higher Na cover- 
age. This is evident from recasting Eq. Q]as 

7 = 70 - ^NaMNa, ( 2 ) 

where 6*N a is the adsorbate coverage. From Eq. [2] it is 
clear that reconstructions with larger #Na are increas- 
ingly favored as //Na increases. For a given value of /^N a 
the reconstruction with the lowest surface energy 7 will 
be realized in an experiment, if it is performed under con- 
ditions of thermodynamic equilibrium. Phase transitions 
can thus occur as ^tNa is changed. 

Thermodynamics places an upper bound on the ad- 
sorbate chemical potential, [i^a, < /^Nai given by the en- 
ergy per atom in the ground-state (body-centered cubic) 
phase of elemental Na. Exceeding this limit in an ex- 
periment would result in precipitation of elemental Na, 
because that phase would then be energetically prefer- 
able to any adsorbed phase. When making the chemical 
potential sufficiently low, by turning the partial pressure 
to a very small value, the bare surface will be the most 
stable phase. The more interesting question is what hap- 
pens in between these two extremes. 

III. RESULTS & DISCUSSION 

A. Monovalent adsorbates 

The calculated DFT surface energies, as a functional 
of chemical potential, are shown in Fig. [2ja) for several 
reconstructions of Si(lll) induced by the monovalent ad- 
sorbate Na. All surface energies are given relative to that 
of the bare Si(lll)-(7x7), which we place at 7=03 The 
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FIG. 2: (Color online) (a) Surface-energy diagram for mono- 
valent Na adsorbates. We compare pure honeycomb chains 
(0=1/3, 1/6 and ML, blue lines) and pure Seiwatz chains 
((9=1/2, 1/4 and ML, green lines). The surface energy of 
Si(lll)-(7x7) (0=0 ML, black line) is also shown, (b) Surface- 
energy diagram for divalent Ca adsorbates. We compare pure 
honeycomb chains (0=1/6, blue line), pure Seiwatz chains 
(0=1/2, green line) and a mixed chain configuration with al- 
ternating honeycomb and Seiwatz chains (0=3/10, red line). 
The vertical dashed line represents the bulk CaSi2 silicide. 



colored lines with non-zero slopes represent the surface 
energies for various Na-induced reconstructions. Results 
for two HCC phases are shown (blue lines), with 9^ a =l/3 
and 1/6 ML corresponding to fully occupied (3x1) and 
half- filled (3x2) channels, respectively. Also shown are 
results for two Seiwatz-chain phases (green lines), with 
#Na=l/2 and 1/4 ML, corresponding to fully occupied 
(2x1) and half-filled (2x2) channels, respectively. Fi- 
nally, results for "empty" HCC and Seiwatz reconstruc- 
tions (i.e. without adsorbates) are shown as flat lines. 

From the energy ordering of these various reconstruc- 
tions, it is evident that over the allowed range of /iNa 
only two phases are thermodynamically stable: the clean 
Si(lll)-(7x7) surface and the (3x1) HCC phase with 
#Na=l/3, in agreement with experiment. We find, but do 
not here show, very similar results for other monovalent 
adsorbates (Li, K), and conclude that for monovalent ad- 
sorbates the only thermodynamically stable phase is the 
HCC with every site in the channel occupied. 



B. Divalent adsorbates 

We turn now to reconstructions induced by divalent 
adsorbates, such as the alkaline earth metals (Mg, Ca, 
Sr, Ba) and rare earth metals (Sm, Eu, Yb). From ex- 
periments we know that at a coverage of 1/6 ML, di- 
valent adsorbates stabilize (3x"2") honeycomb chains. 
At a coverage of 1/2 ML, the (2x1) Seiwatz chains are 
experimentally observed * 21 ! 22 At intermediate coverages, 
divalent adsorbates are known to stabilize mixed chain 
structures with higher periodicity consisting of a com- 
bination of honeycomb chains and Seiwatz chains. The 



simplest combination alternates between honeycomb and 
Seiwatz chains, resulting in (5x1) unit cell, as shown in 
Fig- IHc). The (5x"2") periodicity observed for divalent 
adsorbates may be thought of as being built from two 
(3x1) honeycomb-chain units and two (2x1) Seiwatz- 
chain units. 

We have calculated the surface energies of a variety of 
candidate reconstructions based on the divalent adsor- 
bate Ca, including pure honeycomb chains with cover- 
ages 0= 1/3 and 1/6 ML, pure Seiwatz chains with cov- 
erages 8—1/2 and 1/4 ML, and the mixed configuration 
alternating between honeycomb and Seiwatz chains with 
adsorbate coverages 0=2/10, 3/10, and 4/10 ML. A new 
consideration arises for Ca, because Si and Ca can form a 
variety of stable bulk silicides, such as CaSi2 and Ca2Si. 
To prevent precipitation of these bulk phases, the Ca and 
Si chemical potentials must also satisfy the inequalities 

jUCa + 2/isi < M(CaSi 2 ), (3) 
2Mca + Aisi < Ai(Ca 2 Si), (4) 

where /x(CaSi2) is the energy per formula unit of CaSi2, 
and likewise for Ca2Si. Since /isi is fixed, these con- 
straints have the effect of further lowering the highest 
allowed value of /^ca- 

The resulting surface energies for Ca-induced recon- 
structions are shown in Fig. [2]Jb) . To keep the figure un- 
cluttered we have only plotted the phases that are stable, 
or close to stable, for some allowed value of fica- Two Ca- 
induced phases are thermodynamically stable: the HCC 
reconstruction with 6—1/6 ML (blue) and the Seiwatz- 
chain reconstruction with 9=1/2 ML (green). The mixed 
HCC+Seiwatz configuration with 0=3/10 ML (red line) 
is energetically just above these two phases, but passes so 
close (within 1 meV/A 2 ) to their intersection point that 
its formation cannot be ruled out. Experimentally, the 
mixed configuration is indeed found at coverages between 
those of the two pure phases. 

C. Trivalent adsorbates 

Experiments using trivalent adsorbates (Gd, Dy, Er, 
Ho) reveal a mixed configuration with (5x"2") period- 
icity consisting of alternating honeycomb and Seiwatz 
chains, as shown in Fig. [He). Pure honeycomb or pure 
Seiwatz chain structures are not observed. We have cal- 
culated the surface energies of a number of hypothet- 
ical Gd-induced configurations, including pure honey- 
comb chains with every channel site occupied (0=1/3 
ML), every second site occupied (0=1/6 ML), and ev- 
ery third site occupied (0=1/9 ML), as well as of Seiwatz 
chains with every channel site occupied (0=1/2 ML), and 
every second site occupied (0=1/4 ML). We also consid- 
ered three mixed configurations with 0=2/10, 3/10, and 
4/10 ML. Finally, we also calculated the surface energy 
for the well-studied epitaxial GdSi2 silicide, which con- 
sists of 1 ML of Gd on Si(lll) in the so-called "B-T4" 
structure! 23 ! 24 
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FIG. 3: (Color online) (a) Surface-energy diagram for the 
trivalent adsorbate Gd. We compare pure honeycomb chains 
(0=1/9, blue) and a mixed chain configuration with alternat- 
ing honeycomb and Seiwatz chains (0=2/10, red). The hori- 
zontal line represents Si(lll)-(7x7) (0=0 ML, black), and the 
dashed line represents the GdSi2 silicide film (0=1 ML), (b) 
Surface energies of Gd phases plotted as a function of cover- 
age. Only the points have meaning; the lines are added for 
clarity. 



The resulting surface energies for Gd-induced recon- 
structions are shown in Fig. [3]Ja). As before, we plot 
only the phases that are stable, or nearly so. There are 
two thermo dynamically stable phases, the bare Si(lll)- 
(7x7), and the GdSi2 silicide phase with Oq^—I. None of 
the Gd-chain reconstructions is thcrmodynamically sta- 
ble within DFT. 

To investigate the possibility that the experimentally 
observed phases are kinetically stable, we examine in 
more detail the energetics of all phases with intermediate 
Gd coverage. In particular, we are interested in the ten- 
dency of an adsorbate phase with intermediate coverage 
8cd to separate into two stable phases, one with lower 
coverage Q~^ d and one with higher coverage 0Q d . Let x 
denote the fraction of the total surface area occupied by 
the higher coverage phase. Then the average coverage of 
the two phases is simply x9q d + (1 — x)9' Gd . If x is cho- 
sen such that this average coverage is equal to the cov- 
erage 6cd of the homogeneous phase, then the surfaces 
energies of the homogeneous and inhomogeneous phases 
can be compared directly, without the need for specify- 
ing a chemical potential. To make this comparison simple 
we return to Eq. 2 and now regard the surface energy 
gamma as a function of coverage Ocd for a fixed, arbitrary 
value of chemical potential ficd- Although the surface en- 
ergies for the individual phases depend explicitly on /led, 
the energy difference between the two alternative scenar- 
ios (the homogeneous phase versus the phase-separated 
phase) with the same average coverage does not. For 
display purposes we choose a value for which the two 
endpoint phases have equal surface energies. 

The resulting surface energies, for all the Gd chain 
phases we have considered, are plotted in Fig. (3Jb) rel- 
ative to the endpoint phases, as a function of Gd cover- 
age. For intermediate Gd coverages, the surface energy 




FIG. 4: (Color online) (a) Empty and (b) fiUed-state STM 
image of Si(lll) surface after Gd deposition and subsequent 
annealing. Atomic chains coexist with clean areas with the 
Si(lll)-(7x7) reconstruction separated by triangular GdSi2 
islands. U=±1.8 V, 7=0.18 nA, 45 nmx45 nm. 

is always higher than at the endpoints. Thus, a surface 
prepared with intermediate coverage will, under condi- 
tions of thermodynamic equilibrium, phase separate into 
an appropriate mixture of the two endpoint phases. But 
what about conditions under which thermodynamic equi- 
librium cannot be achieved? 

There are good reasons to consider this scenario. The 
conversion of an intermediate phase to a combination of 
bare (7x7) and the GdSi2 silicide phase requires a con- 
siderable rearrangement of the top several atomic layers; 
this restructuring may be kinetically hindered. On the 
other hand, the structural differences among the various 
Gd-chain phases are relatively minor: for example, all 
are quasi-one dimensional with similar underlying build- 
ing blocks. Thus, the thermodynamically favored end- 
point phases may be kinetically inaccessible, even though 
equilibrium is achieved among the subset of structurally 
similar chain phases with intermediate coverage. 

To analyze the consequences of this hypothesis, we ap- 
ply the Maxwell construction to the phases with inter- 
mediate Gd coverage. This geometrical analysis allows 
one to answer the question of which, if any, intermedi- 
ate phases are stable against phase separation into ap- 
propriate mixtures of other phases. We consider again 
the results in Fig. [3^b) and construct an analogous plot 
(not shown), this time with the (7x7) and GdSi2 end- 
points excluded. Of the remaining phases, only one is 
stable with respect to phase separation into all possible 
other pairs: the (5x2) mixed HCC+Seiwatz phase with 
#Gd=2/10. This conclusion is easy to visualize in Fig. 
|3£b) by simply drawing straight lines between all possible 
pairs of phases; the dotted line shows one such example, 
demonstrating the stability of the #Gd=2/10 phase with 
respect to separation into the two endpoint phases, with 
e Gd =l/9 and G d=l/2. 

These findings are indeed consistent with our exper- 
imental observations. Starting with the deposition of 
2/10 ML Gd onto the substrate held at 680° C, followed 
by short annealing still at 680° C and a cool-down to 
room temperature, one obtains a surface uniformly cov- 
ered with atomic chains. After annealing for a slightly 
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longer period or at slightly higher temperatures, phase 
separation may directly be observed in STM images (see 
Fig. [4} : triangles of GdSi2 silicide form together with re- 
gions of clean Si(lll)-(7x7) coexisting with the chains. 
Further annealing completely transforms the chains into 
silicide islands and areas with Si(lll)-(7x7). 

IV. CONCLUSION 

We used first-principles total-energy calculations and 
scanning tunneling microscopy to study the stability and 
structure of atomic chains of monovalent, divalent, and 
trivalent adsorbates on Si(lll). 

For monovalent adsorbates the theoretical and experi- 
mental results are in excellent agreement, both identify- 
ing the honeycomb-chain channel (HCC) reconstruction 
with adsorbate coverage 9 —1/3 as the only stable chain 
phase. 

For divalent adsorbates three chain phases are found 
experimentally, corresponding to coverages 1/6, 1/2, and 
intermediate values. These findings are corroborated by 
our total-energy calculations, which identify the three 
lowest-energy phases as a half-occupied HCC phase, a 



fully occupied Seiwatz-chain phase, and a simple combi- 
nation of these. 

For trivalent adsorbates, our total-energy calculations 
indicate that the observed combination of HCC and Sei- 
watz phases with adsorbate coverage 2/10 is kinetically 
stable with respect to phase separation into other, ther- 
modynamically stabler phases. 

For all adsorbates, the thermodynamically (or ki- 
netically) stable phases all obey a surprisingly simple 
electron-counting rule proposed earlie r 10 i 15 . 

Acknowledgments 

We thank Franz J. Himpsel for stimulating discus- 
sions. The help of Leslie- Anne Fendt, Samuel Hoffmann, 
Claude Monney and Christoph Walther is gratefully ac- 
knowledged. Skillful technical assistance was provided 
by our workshop and electric engineering team. This 
work was supported by the Fonds National Suisse pour 
la Recherche Scicntifique through Div. II, MaNEP and 
the US Office of Naval Research. Computations were 
performed at the DoD Major Shared Resource Center at 
ASC. 



* Electronic address: corsin.battagliaOunine.ch 
URL: http : //www.unine . ch/phys/spectro 

1 S. C. Erwin, A. A. Baski, and L. J. Whitman, Phys. Rev. 
Lett. 77, 687 (1996). 

2 A. Laracuente, S. C. Erwin, and L. J. Whitman, Phys. 
Rev. Lett. 81, 5177 (1998). 

3 T. An, M. Yoshimura, I. Ono, and K. Ueda, Phys. Rev. B 
61, 3006 (2000). 

4 A. A. Stekolnikov, J. Furthmiiller, and F. Bechstedt, Phys. 
Rev. Lett 93, 136104 (2004). 

5 C. Battaglia, K. Gaal-Nagy, C. Didiot, C. Monney, E. F. 
Schwier, M. G. Garnier, and P. Aebi, to be published 
(2008). 

6 C. Gollazo-Davila, D. Grozea, and L. Marks, Phys. Rev. 
Lett. 80, 1678 (1998). 

7 L. Lottermoser, E. Landmark, D.-M. Smilgies, M. Nielsen, 
R. Feidenhans, G. Falkenberg, R. L. Johnson, M. Gierer, 
A. P. Seitsonen, H. Kleine, et al., Phys. Rev. Lett. 80, 3980 
(1998). 

8 S. C. Erwin and H. H. Weitering, Phys. Rev. Lett. 81, 2296 
(1998). 

9 R. Seiwatz, Surf. Sci. 2, 473 (1964). 

10 C. Battaglia, H. Cercellier, C. Monney, M. G. Garnier, and 
P. Aebi, Europhys. Lett. 77, 36003 (2007). 

11 J. N. Crain, A. Kirakosian, K. N. Altmann, C. Bromberger, 
S. C. Erwin, J. L. McChesney, J.-L. Lin, and F. J. Himpsel, 
Phys. Rev. Lett. 90, 176805 (2003). 

12 J. R. Ahn, P. G. Rang, K. D. Ryang, and H. W. Yeom, 
Phys. Rev. Lett. 95, 196402 (2005). 

13 J. Guo, G. Lee, and E. W. Plummer, Phys. Rev. Lett. 95, 
046102 (2005). 



14 P. C. Snijders, S. Rogge, and H. H. Weitering, Phys. Rev. 
Lett. 96, 076801 (2006). 

15 C. Battaglia, H. Cercellier, C. Monney, L. Despont, M. G. 
Garnier, and P. Aebi, J.Phys.: Conf. Series 100, 052078 
(2008). 

16 We define 1 ML as the atom density of the bulk-terminated 
Si(lll) surface of 7.8 x 10 14 cm" 2 . 

17 G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993). 

18 G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996). 

19 E. Kaxiras, Y. Bar- Yam, J. D. Joannopoulos, and K. C. 
Pandey, Phys. Rev. B 35, 9625 (1987). 

20 The surface energy for the Si(lll)-(7x7) reconstruction 
was determined by calculating the surface energy of the 
Pandey model 2 ^ for Si(lll)-(2x 1) and subtracting the 
published DFT energy difference of 60 meV per (lxl) unit 

ce m. 

21 A. A. Baski, S. C. Erwin, M. S. Turner, K. M. Jones, J. W. 
Dickinson, and J. A. Carlisle, Surf. Sci. 476, 22 (2001). 

22 T. Sekiguchi, F. Shimoksoshi, T. Nagao, and S. Hasegawa, 
Surf. Sci. 493, 148 (2001). 

23 C. Rogero, C. Koitzsch, M. E. Gonzales, P. Aebi, J. Cerda, 
and J. A. Martin-Gago, Phys. Rev. B 69, 045312 (2004). 

24 C. Bonet, I. M. Scott, D. J. Spence, T. J. Wood, T. C. Q. 
Noakes, P. Bailey, and S. P. Tear, Phys. Rev. B 72, 165407 
(2005). 

25 K. Pandey, Phys. Rev. Lett. 74, 1913 (1981). 

26 K. D. Brommer, M. Needles, B. Larson, and J. D. 
Joannopoulos, Phys. Rev. Lett. 68, 1355 (1992). 



